UNCLASSIFIED 

ao  430722 

DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION,  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


NOTICE :  When  government  or  other  drivings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  vlth  a  definitely  related 
government  procurement  operation,  the  U.  3. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  fumlahed,  or  in  any  nay 
supplied  the  said  drawings,  specifications,  or  ether 
data  la  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  In  any  way  be  related 
thereto. 


722 


MEMORANDUM  REPORT 
M64-14-1 


UNITED  STATES 


FRANKFORD  ARSENAL  MEMORANDUM  REPORT  M64-14-1 

Research  and  Development  Group  December  1963 

Pitman-Dunn  Institute  for  Research  OMS  Code  5010.  11.83800.  01 

Philadelphia,  Pa.,  19137  DA  Project  599-01-004 


A  COMPUTER  PROGRAM  FOR  A  KRAMERS -KRONIG 
TRANSFORMATION  OF  THE  OPTICAL  REFLECTIVITY 


Prepared  by: 

R.J.  ESPOSITO 
Physicist 


R.  M.  ROBBINS 
Mathematician 


P.evie-wed  by: 


J.  N.  BROWN 
Mathematician 

syzstJu 

E.  R.  THILO 
Director 

Physics  Research  Laboratory 


Approved  by: 


»-  />-/  >  V. 

A,  L^AylESON 
Depiity  Chief 

Pitman-Dunn  Institute  for  Research 


ABSTRACT 


Relations  involving  the  frequency-dependent  optical  constants 
give  information  on  the  electronic  properties  of  solids.  Power  re¬ 
flectivity  measurements  at  normal  incidence  over  a  wide  frequency 
region  can  be  utilized  in  obtaining  the  optical  constants  by  means  of  a 
general  integral  transformation  known  as  the  Kramer s-Kronig  trans¬ 
formation.  A  computer  program  based  on  this  integral  transformation 
has  been  constructed  in  Fortran  II,  version  9000,  language  for  a 
Remington  Rand  Solid  State  90  computer.  The  program  converts  ex¬ 
perimental  reflectivity  measurements  into  parameters  containing  the 
optical  constants. 

A  brief  presentation  of  the  electronic  quasi-particle  and  collec¬ 
tive  excitation  modes  is  given  in  terms  of  the  optical  constants.  A 
discussion  of  the  approximations  made  in  putting  the  integrals  on  the 
computer  is  also  given.  Two  tests  were  made  to  check  the  accuracy 
of  the  program.  First,  a  test  reflectivity  function,  for  which  the 
Kramers -Kronig  integrals  could  be  directly  evaluated,  was  programmed. 
Then,  published  reflectivity  measurements  on  copper  were  analyzed. 

The  optical  constants  determined  by  the  computer  program  were  in 
good  agreement  with  both  the  test  optical  constants  and  the  previously 
calculated  constants  of  copper. 
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INTRODUCTION 


There  has  been  much  theoretical  speculation  on  the  electronic 
properties  of  alloys.  An  increasing  amount  of  research  is  being  done 
on  their  optical  properties,  elastic  constants,  specific  heat  and  X-ray 
fluorescent  emission  spectra  in  order  to  elucidate  their  band  struc¬ 
tures.  Recently  de  Hass-van  Alphen  oscillations  were  observed  in 
some  IB-IIB  alloys  and  other  alloys.  ^  The  ordered  beta-brass  type 
(IB-IIB)  alloys  are  a  particularly  interesting  family,  for  they  exhibit 
marked  color  changes  as  a  function  of  temperature  in  the  beta-phase 
and  long  electronic  mean  free  paths.  Muldawer^  has  described  the 
color  trend  as  going  from  green  to  yellow,  to  red,  to  gray,  with  in¬ 
creasing  temperature.  For  example,  visual  observations  on  ordered 
beta  CuZii  reveal  that  it  is  yellow  with  a  green  tinge  at  -195*  C,  yellow 
gold  at  28°  C,  and  copper-red  at  300°  C.  These  striking  color  changes 
suggest  that  the  brasses  are  quite  amenable  to  optical  measurements 
as  a  function  of  wavelength  and  temperature.  Also,  their  low  residual 
resistivities,  which  were  observed  by  Rothwarf  and  Muldawer,  3  give 
rise  to  the  possibility  of  making  Fermi  surface  measurements  on  them. 

This  report  will  be  concerned  with  the  construction  of  a  Kramers  - 
Kronig  computer  program  in  order  to  determine  the  electronic  proper¬ 
ties  of  solids  from  the  spectral  reflectivity  at  normal  incidence.  The 
program  has  been  written  by  two  of  the  authors’''  using  Fortran  II,  re¬ 
vision  9000,  for  the  Remington  Rand  Univac  Solid  State  90  computer. 

An  analysis  of  reflectivity  measurements  on  beta  CuZn  is  given  else¬ 
where.  ^ 

The  power  reflectivity  at  normal  incidence  is  not  a  sensitive 
probe  of  the  electronic  behavior  in  a  metal  or  alloy,  but  several  re¬ 
lations  involving  the  optical  constants  give  information  on  the  nature 
of  the  elementary  electronic  excitations.  It  is  well  known  that  elec¬ 
trons  in  a  solid  exhibit  both  individual  and  collective  aspects.  ®  When 
a  fast  charge  particle  passes  through  a  thin  film,  it  can  impart  a 
quantum  of  energy  to  the  valence  electrons.  The  response  of  the  elec¬ 
trons  in  the  valence  band(s)  to  the  high  kinetic  energy  of  the  external 
particle  is  such  that  an  oscillation  in  the  electron  density  occurs.  This 
organized  oscillation  of  the  valence  electrons  as  a  whole  is  an  elemen¬ 
tary  excitation  of  the  electron  gas  and  is  known  as  the  plasma  oscillation. 
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The  free  plasmon  (the  quantum  of  the  plasma  oscillation)  has  an  energy- 
equal  to-fkc'p  where  LC p2  =  4lfNe^/m,  and  N  and  m  are  the  electronic 
density  and  mass,  respectively. 

Plasmons  do  not  influence  the  metallic  behavior  under  ordinary 
circumstances  because  their  excitation  energies  are  at  least  a  few 
electron  volts  (e.  v.  )  greater  than  the  kinetic  energy  of  an  electron  at 
the  Fermi  surface.  In  those  metals  where  the  quasi -particle  core  ex¬ 
citation  energies  are  comparable  to  the  conduction  electron  plasma 
energies,  the  plasmons  are  severely  damped  and  considerably  shifted 
from  their  free  electron  values,  °  For  example,  the  free  electron 
plasmon  in  silver  should  occur  at  approximately  9  e.v.  ,  but  is  actually 
shifted  to  3.  8  e.  v.  ^ 

To  first  order,  it  has  been  shown  that  normally  incident  light 
does  not  excite  a  plasmon  because  the  longitudinal  plasma  modes  do 
not  couple  with  the  transverse  photons.  8*  '  However,  in  the  dielectric 
formulation  of  the  many  body  problem  it  has  been  shown  that,  for  iso¬ 
tropic  materials,  the  longitudinal  dielectric  constant  seen  by  a  fast 
moving  electron  is  equal  to  the  transverse  optical  dielectric  constant 
in  the  random  phase  approximation  at  wavelengths  long  compared  to 
the  average  interelectronic  spacing.  8  Near  the  plasma  frequency,  the 
real  part  of  the  dielectric  constant  approaches  zero,  and  the  imaginary 
part  of  the  inverse  complex  dielectric  constant  (Imag  1/e)  is  a  maximum. 
It  is  just  this  latter  quantity  that  is  proportional  to  the  rate  of  energy 
transfer  from  the  incident  charged  particle  to  the  solid  in  the  electron 
characteristic  energy  loss  experiment.  Hence  the  plasma  effects  may 
be  investigated  in  normal  incident  optical  reflection  and  transmission 
experiments  on  account  of  the  equivalence  of  the  longitudinal  and  trans¬ 
verse  dielectric  constants.  As  the  plasma  frequency  is  approached 
from  low  energies,  a  non-interacting  electron  gas,  due  to  its  mass 
inertia,  is  no  longer  able  to  follow  the  transverse  oscillations  of  the 
normal  incident  electromagnetic  field.  At  this  energy,  a  medium  ex¬ 
hibits  a  characteristic  crossover  from  highly  reflective  to  transparent 
properties  (neglecting  interband  transitions). 

If  a  p-polarized  electromagnetic  wave  (polarized  parallel  to  the 
plane  of  incidence)  is  incident  on  a  metallic  film  at  an  oblique  angle, 
the  surface  plasma  mode  may  be  excited  by  the  electric  field  compo¬ 
nent  normal  to  the  film.  ^  At  the  plasma  frequency  this  interaction 
manifests  itself  by  a  dip  in  the  transmission  or  a  peak  in  the  reflec¬ 
tion.  There  is  no  such  structure  at  the  plasma  energy  for  light  polar¬ 
ized  parallel  to  the  piane  of  the  film. 
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On  the  other  hand,  it  is  well  known  that  photons  excite  single 
electron  intraband  and  interband  transitions.  Structural  changes  in 
the  frequency  dependent  absorption  coefficient  have  been  attributed  to 
direct  interband  transitions  at  critical  points  in  the  joint  density  of 
states  for  the  two  bands  participating  in  the  transition.  ^  These  ob¬ 
served  transitions  at  high  symmetry  points  in  the  Brillouin  zone  serve 
as  a  check  on  the  accuracy  of  theoretical  band  calculations. 

The  transformation  of  the  reflectivity  into  the  optical  constants 
is  a  special  case  of  a  very  general  integral  transformation  known  as 
the  Kramers -Kronig  transformation.  12  When  the  principles  of  line¬ 
arity  and  causality  apply,  the  real  and  imaginary  components  of  a 
frequency  dependent  system  are  interdependent.  That  is  to  say,  if 
the  real  part  of  a  function  is  known  over  a  wide  region  of  the  frequency 
spectrum,  then  the  imaginary  part  may  be  found  by  a  Fourier  integral 
transformation  of  the  real  part,.,  and  vice  versa.  These  dispersion  re¬ 
lations  hold  for  many  linear  systems  whfen  causality  is  not  violated, 
and  have  been  applied  to  such  diverse  fields  as  network  theory,;  the 
quantum  theory  of  scattering,  the  dielectric  formulation  of  the  many 
body  problem,  magnetic  resonance,  etc. 

A  system  is  linear  if  the  response  to  two  or  more  inputs  is  the 
sum  or  superposition  of  each  input;  e.  g.  ,  if 


and  =  c  E2 

then  e  (Ei  +  E2)  =  f  E  =  D  =.  Dj  +  D2 

for  an  isotropic  medium  where  D  is  the  displacement  vector,  €  the  di¬ 
electric  constant,  and  E  the  electric  field. 

The  principle  of  causality  is  concerned  with  a  proper  chrono¬ 
logical  sequencing  of  the  input  and  output  signals.  It  may  be  described 
as  "no  output  can  occur  before  the  input";  "no  scattered  wave  can  ap¬ 
pear  until  the  primary  wave  has  reached  some  part  of  the  scatterer"; 
or  "no  signal  ca.n  be  transmitted  faster  than  the  velocity  of  light  in  a 
homogeneous  material."^  Hence  if 

D  (t)  =  0  at  t  <  tQ  , 
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then  E  (t)  =  0  at  t  <  tQ  . 

Also  the  composition  of  the  system  does  not  change  with  time;  i.  e.  ,  if 

D  (t)  =  €  E  (t) 

— i  — i 

then  D  (t  -  t0)  =  €  E  (t  -  t0). 
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Other  pertinent  relations  are 


c  =  cj  -  ie2  B  (a  -  ik)2 


(5) 


fl 


=  n2 


-  k‘ 


(6) 


(7) 

(8) 

where  e  is  complex  dielectric  constant,  a  the  absorption  coefficient, 
and  X  the  wavelength. 

The  natural  logarithm  of  Equation  1  is 

In  r  =  i  In  R  -  i<f  (9) 


and 


€  2  =  2nk 


a  = 


4ffk 


The  imaginary  part  of  Equation  9  may  now  be  written  as  an  integral 
transform  of  the  measured  reflectivity  by  using  the  Kramers -Kronig 
relation^ 


o 


In  R 


dE 


(10) 


where  E  is  the  energy  of  the  incident  radiation.  Hence,  a  measure¬ 
ment  of  the  reflectivity  over  the  entire  energy  spectrum  enables  one 
to  calculate  the  phase  angle  90  and  the  optical  constants  at  an  arbitrary 
output  energy  E0 . 

In  general,  one  cannot  make  reflectivity  measurements  over  an 
infinite  energy  interval,  for  experimental  difficulties  appear  at  both 
low  and  high  energies.  In  the  soft  X-ray  region  the  reflectivity  be¬ 
comes  quite  small  and  must  be  extrapolated.  The  output  energy  E0 
should  be  far  removed  from  both  the  high  and  low  experimental  limits, 
so  that  the  contributions  beyond  these  energies  to  <£0  are  negligible. 
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Reflectivity  measurements  on  metals  1^  and  semiconductors^  have 
been  carried  out  from  approximately  0.  001  to  25  e.  v. 

In  order  to  treat  the  singularity  of  Equation  10  at  E  =  E0  ,  we 
shall  work  with  the  equivalent  Kramer  a -Kronig  relation^ 


(ID 


where  RQ  is  the  reflectivity  at  E0 .  The  equality  between  Equations 
10  and  11  immediately  follows  since 


Jo 


dE 


and  R0  is  treated  as  a  constant  in  the  integration.  For  the  construc¬ 
tion  of  the  computer  program  it  is  convenient  to  divide  the  contributions 
to  (p  Q  into  different  regions  of  the  energy  spectrum.  Equation  10  will 
be  split  into  four  parts  and  integrated  over  three  intervals,  the  low 
energy  extrapolation  (0  to  a),  the  measured  region  (a  to  b),  and  the 
high  energy  extrapolation  {b  to  ®  ).  Hence, 


where 


<Po  =  <Pl+<P2  +  <P3  +  <P4 


<Pl 


(12) 


(13) 


(T 


2  = 


rr 


In  R0 


+  r°°  dE 

Jb  E2  -  E02 


1 

2rr 


In  R0  In 


(Eo  -  a)  (b  +  Eq) 
(Ec  +  ®)  (b  -  E0) 


(14) 
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The  Integral  <C  i  over  the  measured  region  is  evaluated  numeri¬ 
cally  on  the  computer  by  Simpson's  rule 


fo 

rb  ln  r  ^ 

.  Eo 

ff 

*  a  -  Ee^ 

"  T 

l 3  H 

(Y1  +  4<Y2  +  Y4  +  Y6  + 


.  •  •  +  2{Yj  +  Y5  +  Yy  +  , , ,  Yj,j_2)  +  yn  J  » 


(17) 


where  Y*  is  the  value  of  the  integral  at  the  l4*1  point  from  the  lower 
limit  of  integration  and  AE  is  the  increment  of  E.  Simpson's  rule 
necessitates  that  the  energy  increments  be  equally  spaced  and  that  the 
number  of  points  N  be  an  odd  integer. 


At  the  singularity  point  E  =  E0,  R  =  R0,  the  integrand  becomes 
indeterminant.  By  the  application  of  L'Hospital's  rule,  the  singularity 
may  be  evaluated  in  the  limit  as 


The  derivative 


is  obtained  by  the  approximation  for 


numerical  differentiation  after  smoothing 
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(H)e=Eo  MTJe  [<Rk-2  '  Rr+2)  ‘  8(Rk-i  ’  rk+i}]  (19) 

The  singularity  point  for  a  teat  reflectivity  function  was  found  to  lie  on 
a  smooth  curve  with  the  neighboring  values  of  the  integrand  in  Equation 
17  and  was  included  in  the  Simpson  summation. 

By  inspecting  the  sign  of  the  integrand  in  Equation  13,  one  can  tell 
which  parts  of  the  energy  reflectivity  spectrum  give  a  positive  or  nega- 
tive  contribution  to  the  output  phase  angle.  The  same  is  true  for  Equa- 

2  2 

tion  14.  <p 2  is  positive  when  Ee  ab  and  negative  when  E0  >  ab 
since  In  R0  s  0  as  R  5  1, 

The  reflectivity  in  the  extrapolated  infra-red  interval  (0,  a)  is 
constant  to  a  high  degree  of  approximation;  thus  it  may  be  taken  out 
of  the  integration.  Equation  15  then  reduces  to 

C3.i  (20) 

where  Ra  is  the  long  wavelength  value.  <p  3  gives  a  negative  contribu¬ 
tion  to  the  total  phase  angle. 

The  high  energy  extrapolation  in  many  cases  makes  a  significant 
contribution  to  the  phase  angle.  One  usually  chooses  an  extrapolation 
such  that  the  values  of  the  computed  optical  constants  agree  witli  those 
obtained  from  direct  measurements  in  the  visible  and  the  near  ultra¬ 
violet  frequency  range.  For  simplicity  we  choose  the  reflectivity 
in  the  expression  <£4  to  be  a  constant,  Rc<  Three  different  values  of 
Rc  are  provided  in  this  program.  After  integrating,  Equation  16  be¬ 
comes 

l  /b  +  E„  \ 

p.  =  -  _L  In  Rc  In  ' - -2.)  (21) 

4  2 TT  \  b  -  E0  / 


and  is  always  positive. 
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The  choice  of  a  constant  cut-off  reflectivity  Rc  is  valid  at 
sufficiently  high  energies,  provided  the  range  of  integration  for  0} 
is  large.  In  some  cases  it  may  be  necessary  to  put  extrapolated 
values  in  <fi j  to  get  agreement  with  values  of  the  optical  constants 
measured  by  direct  techniques.  For  this  particular  program  and 
machine,  the  dimensionality  of  the  input  is  500  for  the  <p^  integra¬ 
tion  (actually,  499  input  points  are  available  since  N  must  be  an  odd 
integer).  1®  The  energy  input  points  are  not  entered  as  data  but  are 
generated  by  the  program.  The  starting  point  corresponds  to  E^ 
inputted  as  WA,  with  delta  as  the  increment.  It  has  already  been 
noted  that  energy  increment  must  remain  the  same. 

The  output  points  must  be  appropriately  chosen  input  points 
and  must  be  equally  spaced.  However,  the  output  energy  increment 
need  not  be  the  same  as  the  input  increment.  The  lowest  output 
energy  values  is  £3.  The  other  output,  values  are  determined  by 
chossing  a  desired  increment  (k  =  1 , 2,  3. .  j. 

For  example,  k  -  corresponding  output  points 


1 

E  ,  E  . 

3’.  4’ 

E5  ,,,EN-2 

2 

E3»  *3  » 

E^  , , , 

3 

E3'  E6* 

^  ... 

i 

E3'  E3+i 

*  E3+2i  •** 

The  upper  limit  of  the  output  points  is  not  greater  than  _ 2 *  &n<* 

the  lower  limit  is  not  less  than  Eg.  These  limits  are  due  to  the. 
evaluation  of  the  output  points  near  the  upper  and  lower  limits  of 
the  Simpson  summation.  The  maximum  number  of  output  points  is 
495  and  corresponds  to  the  case  in  which  there  are  499  input  points 
and  k  =  1. 


THE  COMPUTER  PROGRAM 

The  Kramer  s-Kronig  transformation  was  programmed  for 
the  Univac  Solid  State  90  computer,  using  Fortran  II  language, 
version  9000.  The  following  substitutions  were  made  as  a  conven¬ 
ience  of  programming: 
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PHI  =  <fi 


RHO  =  p 
RI  =  R,  R0 

0l  =  <p1 

2  3 

P3=P2+04U)" 

P4=P2  +  4!4  (2) 

^5  =  PZ  +  ^4  (3) 

OMEGA  =••  E,  Eq 

DELTA  =  AE 

WA  =  a  =  Ej 

WB  =  b  =  En 

RC1  , 

RC2 l  =  R.  (Assumed  cutoff  values  at  high  energy) 
RC3,1' 

RA  =  Ra(low  energy) 

CN  =  n 


CK  =  k 


EPSNI  = 

EPSN2  =  e2 
RATI  =  RATIO- 1 


*2 
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RAT  2  =  RATIO-2  = 


WLDO  =  X0  =  L 

The  values  for  R  are  punched  with  nine  columns  per  word 
and  ten  words  to  a  card  in  floating  point.  (If  three  or  four  columns 
make  up  an  input  word,  the  number  of  input  Btorage  locations  would 
still  remain  the  same.)  These  data  values  are  read  by  one  read 
statement  and  must  be  followed  by  a  blank  card. 

Values  for  a,  b,  ®-c;j*  delta,  and  are  punched 

in  floating  points  with  nine  columns  per  word  and  seven  words  on  one 
card.  These., are  read  by  a  single  read  statement.  Integer  values 
for  hi  and  K  are  on  one  card  which  is  read  by  a  single  read  state¬ 
ment. 

The  output  program  prints  the  following  heading:  E#  L,  R, 
PHI,  n,  k,  EPSILON-1,  EPSILON-2<  RATIO-1,  and  RATIO -2  for 
each  of  the  five  values  of  RHO.  After  the  print  out  for  the  first  out¬ 
put  point,  the  paper  is  spaced  and  the  cycle  is  continued  pntil  the 
final  output  is  calculated. 

The  computer  program  is  as  follows. 

C  USS  FORTRAN  It  ***  VERSION  9000 

NO  CAROS 

C  NRAMgRS  KRONta  CALCULATIONS 

DIMENSION  OMESaISOU) .HI! SOU ) • PHI ( 6 ) 1RH0 ( 5I,CN(3),CK<3 ) ' EF5N1 ( 5  I .tPSNJI 51 1 
1KAT 1 ( 5 ) ,  RAT  2(5) 

COMMON  OMEuA.Rl 

Hfc.AU  U.«A>*d*RCi,RCKiHC  J  i  DELTA  >  HA 

READ  14,  N • K 

RE AU  11 >  (KI (J) .  I»1 ,N) 

RCILNRLNIRCI ) 

RC2LN«LN(KC2 ) 
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RCJLHAlN,'RC3) 

RALNALuOU) 

DELIMI./IU.OAUEUTA) 

0£L2«DtLTA/3,0 

GA6>4631854 

SUM2Y«g. 

OMEGA ( 1 )  AWA 
DO  2  IdStN 

2  OMEGA  U)AOM£<iA(I-l)*DEt.TA 
PRINT  19t 
NMiAN-I 

00  7  JdSiNHl *K 
SUM4YA0 • 

W02A0MEGAI J>«*2 
ROLNALNtRI < J) ) 

IG01A1 

IG02A2 

NIA1 

N2H1 

4  00  3  IAN1.N2 

XBOMtGA ( 1 l**2-*02 

aaroln-uhki  (in 

YAA/X 

GO  TO  (33(34) tlGQ2 
33  SUM4YA!sUM4Y*Y 
IG02A2 
GO  TO  3 
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34  SUM2Y»SUN2Y*Y 
1002*1 
3  CONTINUE 

GO  TO  { 8 1 9 • 10 1 •  ISO! 

a  vi#r 

SUM2YAJ • 

M*2 
N2*J-l 
1001*2 
GO  TO  4 

9..N1AJ+1  _  ... 

N2*N 
1001*3 
GO  TO  4 

10  SUM2Y*SUi>12Y-y 

Y*(-l. »/(2.0*Rl(J)*OHEa*tUI !*DELl*l«t I J-2I- (8. 0*H I IJ-1) » ♦ < « .0*N t t J+l > )•*! < 

1  j+2) ) 

joe*j/2 
j 1*2* joe 

1F<J-Jl>  1 7  <  17 1 18 

18  SUM2Y*SUN2Y+Y 
60  TO  19 

17  SUM4Y«SUM*Y*Y 

19  PHI  ( 1 ) A ONE 6 A I J )/3« 1415927*DEL2*t  Yl*4»0*SUH4Y+2 «0*8UH2Y I 
WH» ASOMESA t  J 1 »* A 

*APV*«A+OMEGAI J) 

0*LN(«N*A/»APW 1 
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»BP»*WB*0ME3A( Jl 
»BM»««B-OM£a*( J) 

E«LN(*BP*/»BM»I 

PMnaxtsoLN/apio+et 

PM  I  ( 3 )  *RALN/$.*tN  ( tAP»/»MWA ) 

M»LN(tBM«/tBP») 

PMI  (4)»«ClLN/a*H 
PHI  IS)»RC2LN/q*h 
PHI  (6l«RC3LN/a*H 
KHOu)»Pnm>*pHi<a» 

RH0(2)#RH0(1)*PHI(3) 

RHO<  J)«RH0(2l  +PHI  ( <4  > 

RM0(4)«RH0(2)+PHI(S) 

RHO ( S  > WRHO ( 2 ) +PM 1(6) 

SI3MA«4.*S0RT(Kl (J) ) 

ALUOR 1 • 2396E -Or/OMEQA ( J ) 

WlU01«l2.$bA37l/4lCO 
DO  6  Ifll.5 

2*1.4RIU»-SIGMA*COS(RHO(II  ) 

CMIl'ItiO-RIIJI  )  n 

CK.  (  n#»15HA*SlN(  RHO(  I )  )/l 

EPSNi ( I >  «CN  < I )*»2-CK( I)»*2 

EPSN4(I)»2.*CN(I)*tK(J) 

RATH  Il#ePBN2(i  >/<tPSNl(  I  )**2+EPSN2  ( I )  **2 1 
RAT2 ( I )««LU01»CK( I > 

CONTINUE 

PRINT  12. OMEGA! J) .PUOO.Ri! J). (KHOI I),Cn( I  I .CM  I). EPSNI I  I ) . EPSN2! I ) < RATI ( 1 1 


II* 


1 i AAT2I 1 1 • !#1 »5) 


7  CONTINUE 

11  FORMAT  ( 10E9»3t / ) 

12  FORMAT  (//F9.3,2AE11.3.F9.3.?<axeil.3>/t3lX7<2Aeil.3)/>> 

13  FORMAT ( 7F6»3) 

14  FORMAT (213) 

13  FORMAT  (6*lH£tl2XlHLteAlMR.10*3MPMl«l2XlHN«UXlMR*4X'»MeP»IL0N-lt4X9HEF»lLp 
lN-2*6X7MRATIO-l«6X7HRAT!0-ai 

STOP 
•  ENO 

COMMENTS  ON  THE  ACCURACY  OF  THE  PROGRAM  AND  METHOD 


The  accuracy  of  the  program  was  verified  by  putting  a  known 
test  function  into  the  computer.  The  tabulated  input  reflectivity  was 
evaluated  from  the  expression 

R  =  e"E  (22) 

and  pi  was  integrated  by  Sirpson's  rule  in  the  interval  (a,b),  Suit¬ 
able  approximations  were  made  for  <fl$  and  so  that  no  discontin¬ 
uities  in  the  selected  output  values  appeared.  From  Equations  13 
and  22,  <p^  may  be  directly  integrated 
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The  Univac  computer  values  were  in  agreement  with  hand  calcula¬ 
tions  from  logarithmic  tables  to  the  third  decinv  !  place. 

We  also  ran  Ehrenreich  and  Phillip  s  reflectivity  data  on 
copper^'  ^  up  to  24.4  e.v.  on  the  computer.  We  obtained  excellent 
agreement  with  their  results  for  the  optical  constants  in  the  regions 
from  0  to  10  e.v.  after  a  suitable  value  of  was  chosen  for  the  ^4 
extrapolation. 

The  cut-off  reflectivity  can  drastically  change  the  values  of 
the  optical  constants  determined  by  the  Kramer s-Kronig  program; 
for  as  Rc  goes  to  zero,  <fl 4  approaches  infinity.  Since  the  reflectivity 
cannot  be  measured  at  extremely  high  energies,  .the  Kramer  s-Kronig 
transformation  is  used  in  conjunction  with  values  6f  the  optical  con¬ 
stants  measured  by  direct  methods.  These  include  Drude-type  vari¬ 
ations  employing  polarized  light,  20  the  transmission  interference 
filter  method  of  Schulz,  21,22  aJl(j  a  measure  of  the  unpolari2ed  re¬ 
flectivity  at  four  different  angles.  23  The  average  cut-off  reflectivity 
(more  correctly,  the  average  contribution)  is  found  by  bracketing 
the  known  optical  constants  by  successive  approximations.  Hence,  " 
one  can  give  a  reasonably  good  estimate  of  the  optical  constants  in 
the  regions  beyond  the  energy  interval  covered  by  direct  measure¬ 
ments,  but  well  within  the  reflectivity  limits. 

When  a  material  exhibits  a  sharp  reflection  edge,  the  input 
points  must  be  closely  spaced  in  order  for  the  Simpson  approximation 
to  be  accurate,  The  integration  will  then  extend  throughout  a 
limited  region  of  the  energy  spectrum,  and  the  optical  constants  may 
be  evaluated  only  within  this  region.  The  reflection  edge  in  silver 
drops  from  90  percent  to  less  than  1  percent  between  3,4  and  3.8 
e.  v.  7,  15  if  equal  energy  intervals  are  chosen  every  0.02  e,  v.  , 
the  high  energy  input  would  be  about  10.0  e.v.  Free  electron  plasma 
edges  which  appear  in  the  infrared  for  doped  semiconductors  are 
more  extreme.  The  reflectivity  of  extrinsic  InSB  goes  from  100  per¬ 
cent  to  1  percent  in  an  energy  interval  of  0,  05  e.v.  24  jn  such  cases 
it  is  necessary  to  increase  the  dimensionality  of  the  input  or  to  break 
up  the  1  integration  into  two  different  energy  intervals. 

It  is  interesting  to  note  that  from  our  choice  of  the  phase 
angle  <p  in  Equation  1,  to  be  physically  meaningful, 

0  2  tp  ff.  (24) 
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Otherwise,  the  attenuating  absorption  coefficient  4trk/A  becomes 
negative  and  the  transmitted  light  is  greater  than  the  incident  beam 
at  a  given  wavelength.  This  is  immediately  seen  from  Maxwell's 
equations,  which  give  the  amplitude  of  a  wave  travelling  in  the  x 
direction  through  an  absorbing  medium  as  " 

A  -  ko  exp  Jj^ffi  -  21T1  jj®  -  .  (25) 

The  sign  and/or  magnitude  of  <pQ  may,  therefore,  be  used  as  a  check 
to  detect  whether  a  card  has  been  mispunched  or  the  energy  interval 
chosen  for  the  numerical  integration  is  top  large. 
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